CDm719MENKHOFF 
DIGITAL FILTER COMBINATION FOR INTERPOLATION 



FIELD OF THE INVENTION 

This invention relates to interpolation systems and in particular, to a digital 
filter combination for the discrete-time interpolation of sample values of a sampled 
signal. 

BACKGROUND OF THE INVENTION 

In discrete-time interpolation of sample values of a sampled signal, the 
sample values must represent a band-limited signal or be band-limited by an 
appropriate digital filter. By means of the interpolation, the sampling rate is to be 
changed by an arbitrary numerical ratio, for example, or an intermediate signal value 
is to be calculated while retaining the respective sampling rate. Such' applications 
are in the processing of digital video and audio signals, particularly in digital 
television receivers, video recorders, audio reproduction equipment, or computer- 
based or television-set-based multimedia processing equipment. 

In the simplest case, such an interpolation is a linear interpolation, which 
starts fi-om adjacent given values of the sampling sequence. This simple interpola- 
tion has a disadvantage in that it suppresses signal fi-equencies equal to half the 
sampling rate. In the case of video signals, this reduces the picture sharpness, and 
may even cause clearly visible picture disturbances if such suppression results in 
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interferences in the case of periodic picture contents. 

^^^^IJie- c op e nding pa t en t application , CD I T 12 W , incorporated by reference 
herein, discloses an interpolator for digital signals which circumvents part of these 
disadvantages by being implemented as a digital filter combination which contains, 
in the direction of signal flow, a discrete-time third-order interpolation filter and a 
linear interpolation filter. The discrete-time interpolation filter forms, in addition to 
the existing sample values, an intermediate value exactly in the middle between two 
sampling instants. Through the discrete-time filter, the sampling rate is effectively 
doubled. In the case of the discrete-time filter described, the frequency response in 
the passband for higher fi-equencies can be changed by increasing or decreasing 
particular frequency components by means of a weighting factor. Using a delay 
chain consisting of three delay stages and a logical combination of the differently 
delayed sample values, the individual fi-equency components are formed, which may 
be weighted differently by means of multipliers and are finally combined. By linear 
interpolation, the interpolated value is then calculated fi"om the two closest secondary 
sample values of the new sampling sequence at the desired interpolating instant. To 
perform the linear interpolation, only a subtracter, a multiplier, and an adder are 
needed in the example given, with the mixture ratio of the two sample values being 
controlled by the multiplier and an applied multiplication factor. The multiplication 
factor corresponds to an interpolating instant normalized to the secondary sampling 



rate. 

The filter function of this filter combination, however, does not meet the 
increased requirements placed on digital interpolation. It is therefore an object of the 
invention to provide an improved digital filter for performing an arbitrary temporal 
interpolation of a sampled signal which permits improved picture reproduction, 
particularly with digitized video signals. 
SUMMARY OF THE INVENTION 

This object is attained by providing a digital filter combination for 
mterpolating primary sample values of a sampled signal using an mth-order discrete- 
time filter and a kth-order continuous-time interpolation filter, with m^3 and k ^ 
2; wherein the discrete-time filter forms n secondary sample values firom at least m 
+ 1 primary sample values at equal time intervals, with n ^ m; and the continuous- 
time interpolation filter forms fi-om at least part of the n secondary sample values an 
interpolated value whose temporal position with respect to that of the primary sample 
values is predeterminable by a normalized interpolating instant dp + t^/T, where t^^ 
is the absolute interpolating instant, and T is the period of the primary sampling rate. 
BRIEF DESCRIPTION OF THE DRAWINGS 

The invention and further advantageous features thereof will now be 
explained in more detail with reference to the accompanying drawings, in which: 

Fig. 1 is a schematic block diagram of an embodiment of a digital filter 



combination according to the invention; 

Figs. 2 and 3 are schematic time diagrams of a primary sampling sequence 
and a secondary sampling sequence, respectively; and 

Figs. 4 and 5 are time diagrams illustrating the selection of the sample values 
at different interpolating instants. 
DETAILED DESCRIPTION OF THE INVENTION 

Before embarking on a detailed discussion, the follov^ng should be 
understood. The invention is predicated on the recognition that the object can be 
attained with a digital filter combination consisting of a discrete-time filter for 
calculating fixed intermediate sample values and a continuous-time filter for 
calculating arbitrary intermediate values. The use of such a combination is 
advantageous because the two filter parts can be optimized independently of each 
other. 

If the third-order discrete-time fiher used in the circuit disclosed in copending 
application CDIT-1230 to calculate the fixed intermediate sample values were 
replaced by a fourth- or higher-order filter, this would have the advantage of giving 
a better amplitude response, but the disadvantage that, because of the impulse 
response extending over many sample values, the transient oscillation of the discrete- 
time filter would also extend over many sample values, which would interfere with 
the interpolation. Therefore, the order of the discrete-time filter is not changed. 
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However, instead of the linear interpolation, which corresponds to a first-order 
interpolation, an at least second-order interpolation is used in the continuous-time 
filter. The interpolation polynomial is designed to be easy to implement as a 
monolithic integrated circuit incorporating the discrete-time filter. For the combined 
filter circuit it is irrelevant to what extent some ftinctional units are implemented as 
separate circuit blocks or realized within a processor by means of a program. 

Referring now to Fig. 1 , there is shown an interpolation filter circuit 37 which 
contains a discrete-time filter 1 on the input side and a continuous-time interpolation 
filter 2 on the output side. The two filters are interconnected by a gang switch 3, 
whose control input 3 1 is supplied with a switching signal ps from a controller 4, 
which forms this signal in dependence on a predetermined normalized interpolating 
instant dp. An input 1 1 of discrete-time fiher 1 receives primary sample values sp 
whose bandwidth is limited by the bandwidth of the original signal sg (see Fig. 2) 
and the Nyquist criterion, Input 1 1 is followed by a delay chain 12 consisting of 
three delay stages: 121, 122 and 123, which each provide a delay equal to the period 
T of the primary sample values sp. By different combinations of the primary sample 
values sp tapped from the delay chain 12, the individual frequency components for 
forming secondary sample values ss are formed. The embodiment of the discrete- 
time filter 1 shown in Fig. 1 requires only a single genuine multiplier 1 8 and only 
four delay stages and four adder-subtracters. The circuit configuration with the three- 
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stage delay chain 12 corresponds to a third-order discrete-time filter, which may also 
be implemented in another manner, of course. 

The mput of the first delay stage 121 is connected to the first input of a first 
adder 1 5, and the output of this delay stage is coupled to the first input of a second 
adder 16, whose second input is connected to the input of the third delay stage 123, 
whose output is coupled to the second input of the first adder 15. The outputs of the 
first and second adders 15 and 16 are connected, respectively, to the minuend input 
and the subtrahend input of a first subtracter 17, whose output is coupled through a 
first multiplier 18 to the first input of a third adder 19, whose second input is 
connected to the output of the second adder 16 and whose output is coupled through 
a second multiplier 1 10 to a second output 132 of discrete-time fiher 1 and through 
a fourth delay stage 111, which provides a delay equal to the period T of the primary 
sampling clock, to a fourth output 134 of discrete-time filter 1. The second and 
fourth outputs 132 and 134 provide those interpolated sample values which lie 
exactly between the primary sampling instants, i.e., at the normalized instants dp = 
+ V2 and db = - Vi, respectively, which correspond to the absolute instants t = i,^f2 and 
t = t.5/2, see Fig. 3. 

In addition to the interpolated intermediate values, the secondary sample 
values ss also contain the original sample values sp. Since the discrete-time filter 1 
provides at its outputs at least three adjacent secondary sample values at the same 
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time, the delay chain 12 is tapped between the first and second delay stages 121, 122 
and between the second and third delay stages 122, 123, and the signals available at 
the taps are transferred to a first filter output 131 and a third filter output 133, 
respectively. Thus, four successive secondary sample values ss are available 
simultaneously at the four outputs 131 to 134 of the discrete-time filter 1, and three 
of them are fed to the continuous-time interpolation filter 2 for forming the respective 
interpolated value st. The selection is made by the gang switch 3, whose output 
provides at least a first, second, or third sample value si, s2, s3. The sampling 
instants of the four outputs 131 to 134 correspond in the same order to the four 
absolute instants t., = *-T, t.3/2 = -3T/2, t.2 = -2T, and t.5/2 = -5T/2. The period T* of 
the secondary sampling clock is T* = T/2. 

In a first switch position pi, the gang switch 3 connects the outputs 131 to 
133 to first, second, and third inputs 311, 312, and 313, respectively, of the 
continuous-time interpolation filter 2. In a second switch position p2, connections 
are established between the second to fourth outputs 132 to 134 and the three inputs 
3 1 1 to 3 1 3. The association between the normalized interpolating instants dp + t^/T 
and d = t*u/T* and the respective secondary sample values ss to be tapped is apparent 
firom Figs. 4 and 5. 

The second multiplier 1 10 only caxises the output signal of the third adder 1 19 
to be weighted with the factor 1/2, which may come fi:om a data memory 14, for 



example. If the circuit is implemented in hardware, and binary-coded binary 
numbers are used, the second multiplier 1 10 will be replaced by a simple arithmetic 
shift. The first multiplier 1 8, however, is a genuine multiplier which multiplies by 
a filter coefficient b with which the frequency response in the passband at higher 
frequencies can be slightly accentuated or lowered. The control range of this filter 
coefficient b is relatively small: -3/8 ^ b ^ -1/8. The filter coefficient may be stored 
as a fixed or progranmiable value in a memory 181, for example. In mathematical 
form, the transfer ftmction H(z) of the discrete-time filter 1 between the terminals 1 1 
and 132 can be represented in the complex frequency domain as 
H(z) = b/2 + (z' + X (1 - b) X Vi + (z'^ x b/2) 

The continuous-time interpolation filter 2 in Fig. 1 is designed to perform a second- 
order Lagrange interpolation with the three applied sample values si , s2, and s3. The 
impulse response h(t) in the time domain of such a ftmction is, for example, 
h(t) = (-1/8 + 1 X t/2) x6.i(t) + (-9/8 + 3t-3tV2) x6.,(t-l) + 

+ (45/8-6t+3t2/2)x6.,(t-2) + (.35/8+3t+t2/2) x6.,(t-3) 
where 6.i(t) is a step function whose transition from the 0 state to the 1 state occurs 
at the instant t. This second-order Lagrange polynomial permits the determination 
of an interpolated value at the instant t via three equidistant sample points. The 
interpolated value may vary in time only within a range which extends over half the 
sample-point spacing on the right and the left of the middle sample point; see also the 
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explanations of Figs. 4 and 5. 

The cx)ntinuous-time interpolation filter 2 in Fig. 1 contains the three inputs 
311,312 and 3 1 3 for the simultaneous supply of the first, second, and third sample 
values si, s2, and s3, respectively. The first input 31 1 is connected to the second 
input of a fourth adder 32 and to the minuend input of a second subtracter 33. The 
second input 312 is connected to a third multiplier 34 and to the first input of a fifth 
adder 35. The third input 3 1 3 is coupled to the first input of the fourth adder 32 and 
to the subtrahend input of the second subtracter 33. The output of the third multiplier 
34 is connected to the subtrahend input of a third subtracter 36, whose minuend input 
is fed by the output of the fourth adder 32. The output of the third subtracter 36 is 
coupled to the input of a fourth multiplier 37, v^hose output is connected to the first 
input of a sixth adder 38, whose second input is connected to the output of the second 
subtracter 33. The output of the sixth adder 38 feeds a fifth muhiplier 39, whose 
output is coupled to the second input of the fifth adder 35. The output of the fifth 
adder 35 provides the required interpolated value st. Since the third multiplier 34 
only has to double the second sample value s2, it may be replaced by an arithmetic 
shift. The fourth and fifth multipliers 37, 39, however, perform genuine multiplica- 
tions, namely by the interpolating instant d normalized to the period T* and by one 
half of the value of this instant, d/2. This half value d/2 can be formed either by 
means of an arithmetic shift or by means of a fiuther multiplier 310, whose two 
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inputs are supplied with the values d and V^, the latter from a data memory 311. 

The controller 4 contains a comparator for forming the switching signal ps 
for the gang switch 3. For this, either the absolute interpolating instant t^^ or the 
normalized interpolating instant dp is evaluated. Furtheraiore, the controller 4 forms 
the interpolating instant d referred to the two closest secondary sample values ss, sp* 
and the associated period T* from the applied interpolating instants t^^ or dp, with d 
= t*i/r*. The two normalized interpolating instants dp and d are henceforth called 
the primary and secondary interpolating instants dp and d, respectively. The 
relationship between the switching signal ps and the primary and secondary 
interpolating instants dp and d will become apparent from the description of Figs. 4 
and 5. While, the meaningfiil interpolation, the range for the primary interpolating 
instant dp extends from -1 to +1 at best, the range for the secondary interpolating 
instant d, which extends from -0.5 to +0.5, can only cover subranges thereof which 
adjoin each other without overlap. According to the first and second switch positions 
pi, p2, the resulting interpolation ranges pi, p2 are given in Fig. 5 beside the time 
axis for the primary interpolating instants dp. 

If only negative or positive time ranges are permitted for the primary 
interpolating instant dp, they will extend between 0^ dp^ 1 or -1 £ p^ 0. In that case, 
too, any intermediate time value between two primary sampling instants to, t.j, t.2, ... 
can be implemented. The respective ranges for the switching signals ps which 
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belong together are apparent from the switch positions pi, p2 given in Fig. 5. 

For any interpolating instant t^n within a period T, however, interpolation is 
also possible if a symmetrical time interval -0,5 ^ dp ^ +0.5 or an unsymmetrical 
time interval -0.25 ^ dp ^ +0.75 is selected for the interpolating instant dp. In the 
first case, three different switch position ranges pi, p2, pi are necessary; in the 
second case, only the two switch-position ranges pi, p2 are required. Each switch 
position is assigned different secondary sample values ss, namely three groups of 
three in the first case and two groups of three in the second case. In the case of a 
limited interpolation interval, a smgle group of three may be suflBcient. In that case, 
the gang switch 3 can be dispensed with. 

The controller 4 forms the secondary interpolating instant d from the primary 
interpolating instant dp either by obtaining the latter from a stored table or by a 
simple modulo calculation d = dp +/-j, with j = 1 , 2, 3, ... having to be selected so that 
the range -0.5 ^ d ^ + 0.5 will not be exceeded. In defining all range limits, 
ambiguities must be avoided, with the conditions of the respective numerical system 
having to be taken into account. In the two's complement system, the positive range 
of values, for example, excludes the power of two, 2", as attainable limit values. 

The schematic time diagram of Fig. 2 shows the analog signal sg, which is 
sampled at instants t.3 to \ . The resulting primary sample values sp represent a 
digitized data sequence. The period T of the primary sampling rate is constant. 
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Therefore, the time axis can also be normalized to the period T, whereby normalized 
time values dp are formed, to which correspond the subscript values of the times t in 
Fig. 2, The formation of normalized time values dp is of advantage for the formation 
of the required interpolation value st. Fig. 2 shows temporally equidistant 
interpolated values stl, st2, st3, which may be assigned to a new sampling sequence, 
for example. The associated period Ts may be in an arbitrary numerical ratio to the 
original period T. 

The time diagram of Fig. 3 shows a sequence of secondary sample values ss, 
which may be formed by the discrete-time filter 1 of Fig, 1, for example. Compared 
with the primary sampling sequence, the period T* of the secondary sampling 
sequence is halved, i.e., T* = T/2. The secondary sampling sequence contains the 
original sample values sp, which are denoted in Figs. 3 to 5 by sp*. 

With respect to the signal sg, the interpolating instant t*^ in Fig. 3 is identical 
to the interpolating instant ti„ of Fig. 2. The is to indicate that in the case of the 
secondary sampling sequence, the interpolating instant has to be referred to the 
respective closest secondary sample value ss, which in Fig. 3 is identical with the 
primary sample value sp*. 

For the interpolation in the continuous-time interpolation filter 2 of Fig. 1, 
groups of three adjacent secondary sample values ss are processed, the selection 
being made such that the interpolating instant t*i„ is as close to a secondary sampling 
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instant as possible. Examples of such a secondary group of three sample values ss 
are shown in Figs. 4 and 5, with the permissible interpolation range t^j^ in Fig. 4 
being symmetrical about the secondary sample value ss at the instant t.3/2, and that in 
Fig, 5 being symmetrical about the primary sample value sp* at the instant t.j. 

In the circuit example of Fig. 1, the selection of the sampling instants 
according to Fig. 4 or Fig. 5 is made by the gang switch 3 in conjunction with the 
controller 4. The temporal limitation of the respective mterpolation range t+j^ 
follows fix>m the interpolation function of the continuous-time filter 2. Interpolated 
values st which lie outside the optimum interpolation interval - 1/2 ^ d ^ + V2 have 
greater deviations from the actual signal waveform sg than those interpolated values 
st which lie within that interval. It is therefore advantageous to use only those three 
sample values ss of the secondary sampling sequence which are closest to the 
interpolating instant t^^ . 

The time diagram of Fig. 4 corresponds to switch position pi of gang switch 
3, and the time diagram of Fig. 5 corresponds to switch position p2. Shown below 
the time axis of Fig. 5 is a further time axis which corresponds to the primary 
interpolating instant dp. For the practical application of the interpolation circuit 
according to the invention, it may be desirable for the interpolating instant t^^ to 
extend over two periods T, namely from a predetermined sampling instant to the 
preceding and subsequent sampling instants. This corresponds to the range of values 
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from -1 to +1 for the primary interpolating instant dp. Above the associated time 
axis dp in Fig. 5, the associated switch positions pi, p2 are given. The associated 
secondary sample values ss are only shown for the time ranges of Fig. 4 and Fig. 5. 
Outside these time ranges, other triples of sample values ss must be used for the 
continuous-time interpolation in filter 2. 
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